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Hydrodynamic  Ram  (HRAM)  is  a  phenomenon  that  occurs  when  a  high-kinetic  energy  object  penetrates 
a  fluid-filled  container.  The  projectile  transfers  its  momentum  and  kinetic  energy  through  the  fluid  to  the 
surrounding  structure,  increasing  the  risk  of  catastrophic  failure  and  excessive  structural  damage.  This  is 
of  particular  concern  in  the  design  of  wing  fuel  tanks  for  aircraft  since  it  has  been  identified  as  one  of 
the  important  factors  in  aircraft  vulnerability.  Usually  the  HRAM  phenomenon  is  analyzed  considering 
completely  filled  tanks,  but  its  effect  on  partially  filled  containers  should  also  be  taken  into  account  due 
to  the  fact  that  tanks  use  to  be  impacted  under  these  conditions.  In  the  present  paper,  the  commercial 
finite  element  code  LS-DYNA  has  been  used  to  simulate  an  HRAM  event  created  by  a  steel  spherical 
projectile  impacting  a  partially  water-filled  aluminium  square  tube.  The  ALE  formulation  is  employed 
to  reproduce  the  event.  Experimental  tests  which  indicate  the  pressure  at  different  points  of  the  fluid, 
displacement  of  the  walls  and  cavity  evolution  for  different  impact  velocities,  are  compared  with  the 
numerical  results  in  order  to  assess  the  validity  and  accuracy  of  the  ALE  technique  in  reproducing  such  a 
complex  phenomenon. 

©  201 1  Elsevier  Masson  SAS.  All  rights  reserved. 


1.  Introduction 

The  process  by  which  a  high-speed  projectile  penetrates  a  fluid- 
filled  tank  and  transfers  kinetic  energy  to  the  surrounding  walls  is 
known  as  Hydrodynamic  Ram  (HRAM).  The  HRAM  effect  in  fuel 
tanks  is  identified  as  one  of  the  important  factors  in  aircraft  vul¬ 
nerability  since  the  fuel  tanks  represent  the  largest  exposed  area 
of  all  the  vulnerable  components.  HRAM  is  especially  dangerous 
for  aircraft  with  extremely  lightweight  designs,  since  the  struc¬ 
tural  resistance  of  their  integral  fuel  tanks  cannot  be  improved 
by  strengthening  the  airframe  because  it  would  counteract  the  re¬ 
quirements  of  a  lightweight  design. 

Vulnerability  to  HRAM  has  been  usually  related  to  military  air¬ 
craft,  but  commercial  airplanes  are  not  exempt  of  its  effect.  In  1990 
the  Federal  Aviation  Administration  (FAA)  established  the  Aircraft 
Catastrophic  Failure  Prevention  Research  Program,  in  which  the 
analysis  of  the  effects  of  an  uncontained  turbine  engine  fragment 
penetrating  aircraft  fuel  tanks  [36]  was  one  of  the  research  areas. 
In  2000,  when  a  Concorde  crashed  after  takeoff  from  Charles  de 
Gaulle  airport  (France),  the  final  investigation  report  revealed  that 
the  HRAM  had  played  a  significant  role  in  the  aircraft  failure. 

Hydrodynamic  Ram  can  be  divided  into  four  principal  stages: 
shock,  drag,  cavitation  and  exit  (Fig.  1).  Each  stage  contributes 
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to  structural  damage  in  a  different  way  and  to  a  different  ex¬ 
tent.  When  the  projectile  penetrates  the  wall  of  the  fluid-filled 
structure,  it  creates  damage  in  the  vicinity  of  the  impact  point; 
in  addition  the  impact  energy  is  transferred  to  the  fluid,  gener¬ 
ating  a  high-pressure  hemispherical  shock  wave.  During  the  drag 
phase,  the  projectile  travels  through  the  fluid,  while  its  kinetic  en¬ 
ergy  is  partially  transformed  into  fluid  motion  as  the  projectile  is 
slowed  by  fluid  drag  forces.  The  displacement  of  the  fluid  from  the 
projectile  path  generates  a  radial  pressure  field.  In  contrast  to  the 
pressure  field  that  develops  during  the  shock  phase,  the  fluid  is  ac¬ 
celerated  gradually  instead  of  impulsively.  This  causes  less  intense 
peak  pressures  but  of  greater  temporal  extent.  The  displacement 
of  fluid  during  the  drag  stage  forms  a  cavity  behind  the  projectile. 
The  subsequent  expansion  and  collapse  (oscillations)  of  the  cavity 
is  known  as  the  cavitation  stage.  The  oscillations  of  the  cavity  can 
cause  significant  pressure  pulses.  The  final  stage  occurs  when  the 
projectile  exits  the  container.  In  contrast  to  the  perforation  of  the 
front  wall,  the  exit  of  the  projectile  occurs  through  a  pre-stressed 
wall,  caused  by  the  initial  shock  stage  and  the  subsequent  loading 
by  the  fluid. 

Simulation  of  HRAM  events  has  been  attempted,  with  more  or 
less  success,  for  over  30  years.  The  first  methods  employed  to  sim¬ 
ulate  HRAM  were  based  on  the  use  of  the  Piston  Theory  for  the 
fluid-structure  interaction.  This  theory  assumes  the  normal  reflec¬ 
tion  of  pressure  waves  when  reaching  the  walls  of  the  structure, 
resulting  in  a  one-dimensional  response  mechanism.  The  Piston 
Theory  was  incorporated  into  two  structural  analysis  codes  in  or- 
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Fig.  1.  Phases  of  Hydrodynamic  Ram. 


der  to  simulate  the  HRAM  phenomenon  [4,2].  The  comparison  of 
the  numerical  results,  provided  by  these  two  codes,  with  exper¬ 
imental  data  [3,11,6]  showed  a  complete  lack  of  agreement  and 
proved  the  failure  of  the  Piston  Theory  for  Hydrodynamic  Ram 
pressure  loading. 

Once  the  Piston  Theory  was  discarded  to  predict  HRAM  fluid- 
structure  interaction,  Lundstrom  [28]  proposed  the  Variable  Image 
Method.  Lundstrom  described  the  flow  field  in  terms  of  a  potential 
function  0  which  satisfied  the  wave  equation.  The  model  attempts 
to  approximate  the  effect  of  the  projectile  and  cavity  by  a  line 
of  sources  distributed  along  the  trajectory  behind  the  projectile. 
This  allows  determining  the  incident  pressure  wave  magnitude  by 
means  of  the  Bernoulli’s  equation  and  then  uses  it  to  calculate  the 
pressure  wave  reflected  from  the  structures’  walls  [29].  In  general, 
this  methodology  has  not  been  particularly  successful  and  the  ge¬ 
ometries  that  can  be  analyzed  are  limited.  However,  the  Variable 
Image  Method  solved  some  of  the  planar  assumptions  in  the  Piston 
Theory  model,  and  provided  a  more  realistic  coupling  between  the 
fluid  and  the  structure.  Other  researchers  [47]  developed  a  plastic 
deformation  model  to  predict  the  structural  response  of  fuel  cell 
walls  due  to  HRAM  in  contrast  to  other  works  in  which  only  the 
elastic  regime  was  considered.  The  model  requires  the  knowledge 
of  the  total  impulse  imparted  by  the  fluid  to  the  fuel  cell. 

Other  codes  such  as  HRSR  (Hydraulic  Ram  Structural  Response) 
[15],  ERAM  or  EHRSR  were  developed  [10],  but  all  of  them  showed 
their  limitations  and  lack  of  accuracy  on  predicting  the  conse¬ 
quences  of  an  HRAM  event  since  none  of  them  fully  coupled 
the  mechanisms  of  fluid-structure  interaction,  nor  did  they  allow 
for  a  complex,  engineered  structure.  The  complicated  physics  and 
mechanics  of  HRAM  phenomena  were  not  satisfactory  solved  un¬ 
til  higher-order  numerical  algorithms  were  incorporated  into  the 
codes  in  the  late  1980’s. 

Coupled  Euler-Lagrange  methods  have  been  under  development 
since  the  early-to-mid  1990’s.  They  combine  the  desirable  charac¬ 
teristics  of  Lagrangian  and  Eulerian  formulations.  These  methods 
are  used  in  multiple  industries  for  a  wide  variety  of  analysis  in 
which  fluids  interact  with  structures  or  when  high  distortions  may 
appear  [36,1,37,16,38,9,7,27,21],  including  airbag  and  tire-water 
dynamics  in  the  automobile  field  [30,22],  the  impact  of  bird  strikes 
on  aircraft  [14,17],  and  the  effects  of  sloshing  on  ships  [50].  The 
Arbitrary  Lagrangian  Eulerian  (ALE)  technique  is  widely  used  in 
those  kind  of  problems.  It  can  be  viewed  as  a  hybrid  between  the 
Lagrangian  and  Eulerian  methods  which  combines  the  best  char¬ 
acteristics  of  both.  It  avoids  the  mesh  distortion  problem  of  the 
Lagrangian  formulation  and  allows  to  precisely  define  the  contours 
unlike  in  the  Eulerian  formulation. 


In  the  last  few  years,  there  have  been  new  advances  in  de¬ 
velopment  and  use  of  computational  methods  for  fluid-structure 
interactions  due  to  the  interest  of  reach  more  effective  computa¬ 
tional  techniques  [5,42,46,20,41,25]  and  solve  more  difficult  prob¬ 
lems  motivated  by  different  industries,  such  as  aeronautics,  naval 
or  more  recently  biomedical  sciences. 

As  an  example  of  the  increasing  interest  on  solving  industrial 
fluid-structure  problems,  it  is  worth  to  mention  the  recent  works 
of  Petitpas  et  al.  [32]  and  Lecysyn  et  al.  [23,24]  in  which  a  ballistic 
impact  on  an  industrial  tank,  filled  with  a  toxic  fluid  and  made 
of  steel,  is  studied.  The  authors  propose  an  analytical  model  to 
reproduce  the  behaviour  of  the  projectile  on  the  fluid  and  study  its 
influence  on  the  toxic  liquid  ejection  and  the  droplets  generated. 
Interesting  results  are  obtained  in  those  works. 

The  modelling  of  coupled  problems  of  fluid-structure  interac¬ 
tion  such  as  HRAM,  has  been  proven  to  be  a  complicated  task  and 
is  still  quite  challenging.  This  is  of  great  importance  since  ana¬ 
lytical  solutions  can  provide  only  a  limited  understanding  of  the 
nature  of  the  behaviour.  In  addition,  the  HRAM  phenomenon  has 
usually  been  studied  in  completely  filled  containers  without  taking 
into  account  its  effect  on  the  structure  on  a  partially  filled  tank. 
Nevertheless,  it  is  more  probable  that  an  aircraft  fuel  tank  could 
be  impacted  when  it  is  not  completely  filled  because  they  are  fully 
filled  for  a  short  period  of  time. 

The  authors  of  this  work  have  previously  observed  the  different 
behaviour  between  completely  and  partially  filled  tanks  impacted 
by  a  projectile  [43].  The  partially  filled  tubes  show  some  particu¬ 
larities  due  to  the  presence  of  air  above  the  fluid  inside  the  tank. 
This  fact  makes  that  the  cavity  inside  the  tank  raises  a  layer  of 
fluid  that  impacts  the  wall.  This  impact  affects  the  structure  in 
a  different  way  than  in  a  completely  filled  tank.  In  addition,  the 
distribution  of  the  fluid  inside  the  tank  influences  the  wall  defor¬ 
mations.  The  authors  also  showed  [44]  the  capabilities  of  different 
numerical  techniques  for  the  simulation  of  the  HRAM  phenomenon 
for  completely  fluid  filled  tanks. 

In  the  present  paper,  numerical  simulations  of  a  partially  water- 
filled  aluminium  square  tube  subjected  to  impact  by  steel  spherical 
projectiles  at  different  velocities  are  shown.  The  simulations  are 
performed  with  the  software  LS-DYNA  applying  the  Arbitrary  La¬ 
grangian  Eulerian  (ALE)  technique  for  the  fluid.  Experimental  tests 
providing  the  pressure  in  different  points  of  the  fluid,  deformation 
of  the  walls  and  cavity  evolution  for  different  impact  velocities  are 
compared  with  the  numerical  results.  The  ALE  technique  was  al¬ 
ready  proved  to  be  valid  on  simulating  the  HRAM  phenomenon  on 
fluid  filled  tanks  [44],  however  it  would  be  necessary  to  confirm 
that  it  also  reproduces  the  particular  characteristics  observed  in 
the  partially  filled  cases.  Under  these  conditions  two  different  flu¬ 
ids  are  present  (water  and  air),  and  hence  the  water  suffers  much 
great  displacements. 

2.  General  description  of  the  ALE  technique 

Lagrangian  methods  are  mostly  used  for  problems  in  which  de¬ 
formation  is  not  extremely  large.  The  mesh  moves  in  space  with 
the  material,  so  that  tracking  surfaces  and  applying  boundary  con¬ 
ditions  can  be  done  easily.  If  the  solids  suffer  large  deformations 
and  mesh  distortion  is  high,  the  calculation  will  stop  and  the  nodes 
of  the  mesh  will  need  to  be  repositioned  in  order  to  continue  the 
calculation.  This  process  is  necessary  every  time  the  mesh  becomes 
too  distorted  for  the  calculation  to  continue.  Therefore,  obtaining 
solutions  to  the  fluid  dynamics  aspect  of  a  high  velocity  impact  in 
a  fluid-structure  interaction  event  (as  HRAM  phenomenon)  is  not 
reasonable  using  a  purely  Lagrangian  method  due  to  the  high  de¬ 
formation  that  the  fluid  mesh  suffers. 

On  the  other  hand,  purely  Eulerian  descriptions  are  typically 
used  for  fluid  calculations.  The  mesh  remains  fixed  in  space  while 
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the  material  passes  through  it.  Surfaces  and  boundary  conditions 
are  difficult  to  track,  however  the  mesh  distortion  is  not  a  problem 
because  the  mesh  is  not  distorted.  Therefore,  obtaining  solutions 
to  the  structural  dynamics  aspects  of  a  high  velocity  impact  in 
a  fluid-structure  interaction  event  is  not  feasible  using  a  purely 
Eulerian  method  due  to  the  deformation  that  the  structural  mesh 
suffers.  The  ALE  formulation  contains  both  pure  Lagrangian  and 
pure  Eulerian  formulations  being  able  to  address  both  the  fluid  and 
structural  dynamics  of  an  event. 

The  ALE  approach  is  based  on  the  arbitrary  movement  of  a 
reference  domain  which,  additionally  to  the  common  material  do¬ 
main  and  spatial  domain,  is  introduced  as  a  third  domain  [39].  The 
problem  is  formulated  in  this  reference  domain,  which  later  will 
correspond  to  the  finite  element  mesh.  The  arbitrary  movement  of 
the  reference  frame  enables  to  deal  with  moving  boundaries,  free 
surfaces  and  large  deformations,  and  interface  contact  problems. 
In  the  ALE  description  of  LS-DYNA  [26],  an  arbitrary  referential  co¬ 
ordinate  is  introduced  in  addition  to  the  Lagrangian  and  Eulerian 
coordinates.  The  ALE  equations  are  derived  by  substituting  the  re¬ 
lationship  between  the  material  time  derivative  and  the  reference 
configuration  time  derivative: 


9/(Xi,t) 

at 


dfQChO 

at 


+  (Vi  -  Ui) 


3/(Xj,t) 

dXi 


9/  (xi»  t) 

at 


+  cot 


9/ (xi.  t) 

dXi 


(i) 


where  X;  is  the  Lagrangian  coordinate,  i  the  referential  coordinate, 
Xi  the  Eulerian  coordinate,  and  v\  and  tq  are  the  material  and  the 
mesh  velocities,  respectively.  In  order  to  simplify  the  equations,  the 
convective  velocity  term  coi  =  u;  -  tq  is  introduced. 

Thus,  the  governing  equations  for  the  ALE  formulation  are  given 
by: 


•  The  conservation  of  mass  equation. 

dp  dv  dp 

-  =  —p - COi  — 

dt  dXi  dXi 

•  The  conservation  of  momentum  equation. 

dVi  dVi 

plF  +  p‘,‘Wl=ai>p  +  pb‘ 

•  The  conservation  of  total  energy  equation. 

dE  dE 

p~di  +  pc0idxj  =  aiiViJ  4  bjVj 


(2) 

(3) 

(4) 


where  p  is  density,  crq  the  stress  tensor,  b[  is  body  force  and  E  is 
energy. 

Note  that  the  Lagrangian  and  Eulerian  formulations  for  the  con¬ 
servation  equations  are  obtained  setting  the  mesh  velocity  equal  to 
the  particle  velocity,  u  =  v,  or  setting  the  mesh  velocity  equal  to 
zero,  u  =  0,  respectively. 

In  the  LS-DYNA  code,  the  ALE  algorithm  is  implemented  using 
the  Split  Operator  Technique  [13]  which  is  a  two  step  process.  The 
ALE  method  begins  working  as  a  Lagrangian  finite  element  code  in 
which  the  motions  of  the  fluid  and  the  structure  are  determined 
and  the  finite  element  grid  is  deformed  following  the  material.  As 
long  as  the  mesh  distortion  is  acceptable,  the  Lagrangian  calcula¬ 
tion  continues.  If  it  becomes  highly  distorted,  the  second  step  (the 
advection  step)  is  performed.  In  the  ALE  methodology  the  distorted 
mesh  can  be  partially  restored  to  its  original  shape  based  on  pre¬ 
defined  criteria  for  element  deformation.  Following  a  prescribed 
measure  of  permissible  distortion  in  the  element,  the  element 
shape  is  changed;  and  mass,  momentum  and  energy  are  fluxed 
across  the  old  element  boundaries  to  calculate  their  new  values  for 


the  new  element  shape.  In  general,  the  ALE  methodology  permits 
flow  of  material  across  element  boundaries  preserving  a  balance 
for  the  physical  velocity  between  grid  motion  and  flow  [10]. 

3.  Experimental  set-up 

In  order  to  achieve  an  appropriate  and  wide  validation  of  the 
numerical  modelling,  it  is  necessary  to  have  enough  experimental 
data.  As  there  is  not  much  available  data  on  pressure,  tank  walls 
deformation  or  cavity  evolution  of  the  HRAM  phenomenon  in  the 
literature,  experimental  tests  were  made  by  the  same  workers  and 
their  results  were  presented  in  a  previous  paper  [43].  The  tests 
consisted  of  high  velocity  impacts  of  steel  spheres  into  tanks  filled 
at  different  levels.  Water  was  selected  as  fluid  for  safety  reasons 
and  because  of  its  similar  density  compared  to  that  of  kerosene. 
The  sketch  of  the  experimental  device  used  for  impact  tests  is 
shown  in  Fig.  2. 

The  test  boxes  consisted  of  6063-T5  square  aluminium  tubes 
750  mm  long,  150  mm  wide  and  2.5  mm  thick.  The  specimens 
were  closed  with  two  PMMA  windows  30  mm  thick,  fixed  to  the 
specimen  with  four  steel  bars;  these  transparent  panels  allow  the 
recording  of  the  impact  process  by  means  of  a  Photron  Ultima 
APX-RS  digital  high-speed  camera.  A  similar  set-up  was  proposed 
by  Nishida  et  al.  [31].  To  obtain  optimal  images  of  the  penetra¬ 
tion  process,  it  is  necessary  an  appropriate  lighting,  which  was 
provided  by  an  Arrisun  12  Plus  lamphead  with  a  1200  W  Hydrar¬ 
gyrum  Medium-arc  Iodide  (HMI)  lamp.  The  contact  points  between 
PMMA  windows  and  the  specimen  were  sealed  with  silicone  in  or¬ 
der  to  avoid  fluid  leakage  [40]. 

Two  pressure  transducers,  PCB  138A06  [40],  were  located  at 
different  points  to  obtain  pressure  data  in  the  fluid.  One  of  the 
transducers  was  placed  near  the  impact  point,  PTn,  and  the  other 
was  placed  far  from  the  impact  point,  PTf.  Two  holes  were  made 
on  the  lower  wall  of  the  specimen  to  place  both  sensors  inside  the 
water.  Their  positions  are  shown  in  Fig.  3.  A  Dewetron  DEWE-800 
data-acquisition  device  was  used  to  record  the  signals.  This  system 
can  record  at  a  sampling  rate  of  1  ps  and  synchronize  the  data 
with  the  video  recording. 

A  one-stage  light  gas  gun  with  a  4.7  litre  chamber,  which  stores 
gas  at  a  maximum  pressure  of  300  bar,  was  used  with  helium.  The 
length  of  the  barrel  was  4.5  m  and  its  calibre  25  mm.  The  gun 
was  aimed  at  an  armoured  steel  specimen  chamber  box  1  x  1  x 
1  m3  inside  of  which  the  specimen  was  placed  during  the  impact 
test.  The  chamber  had  a  small  circular  window  in  the  front  for 
the  projectile  to  pass  through,  and  two  large  lateral  windows  to 
illuminate  the  specimen  and  capture  the  video  sequence  of  the 
impact,  Fig.  2. 

The  projectile  launched  against  the  box  was  a  steel  sphere  with 
a  diameter  of  12.5  mm  and  a  mass  of  8  g.  Two  impact  velocities 
were  chosen:  600  and  900  m/s.  The  test  boxes  were  filled  with 
water  at  different  volume  fractions:  60,  75  and  100%.  In  this  paper, 
only  conditions  corresponding  to  partially  filled  tanks  (60  and  75%) 
will  be  considered  for  the  simulations. 

4.  Numerical  analysis 

The  numerical  model  was  developed  with  the  commercial  fi¬ 
nite  element  code  LS-DYNA  v.971  [26].  This  software  is  particularly 
suitable  for  nonlinear  dynamic  problems,  such  as  impacts  or  explo¬ 
sions.  It  also  allows  the  employment  of  different  techniques  such 
as  ALE  or  SPH  to  solve  fluid  and  fluid-structure  problems.  A  com¬ 
parison  and  validation  of  both  techniques  to  simulate  the  HRAM 
phenomenon  in  a  completely  filled  tube  was  carried  out  by  the 
same  workers  [44].  In  order  to  reproduce  the  HRAM  phenomenon 
in  a  partially  filled  tank,  the  ALE  technique  was  adopted  to  model 
the  fluid  inside  the  tank. 
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Fig.  2.  Sketch  of  the  experimental  device  used  for  the  impact  tests. 


4.1.  Box  and  projectile  FE  model 

The  symmetry  of  the  problem  under  consideration  allowed 
modelling  half  of  the  whole  geometry  (Fig.  4,  left).  Since  the  na¬ 
ture  of  this  simulation  demands  a  very  high  mesh  density,  such  a 
reduction  in  the  model  size  is  very  desirable.  If  the  tube  would  be 
completely  filled  with  water,  the  symmetry  would  allow  modelling 
only  a  quarter  of  the  whole.  The  box  has  been  divided  in  three 
parts,  the  walls  impacted  by  the  projectile  (entry  and  exit  walls), 
the  side  walls  and  the  PMMA  window. 

The  impacted  walls  and  the  PMMA  window  were  discretized  by 
means  of  eight-node  solid  hexahedron  Lagrangian  elements  with 
reduced  integration.  A  refined  mesh,  corresponding  to  the  impact 
zone,  and  a  progressively  coarser  mesh  as  the  distance  to  the  hit 
point  grows,  can  be  observed  in  Fig.  4,  right.  The  impacted  walls 
present  5  elements  through  the  thickness  and  elements  of  1  mm  in 
size  in  the  other  two  directions  near  the  impacted  zone.  Based  on 
previous  simulations,  the  mesh  size  was  considered  appropriate  to 
reproduce  the  behaviour  of  the  solids  in  the  impacted  zone.  Four- 
node  Belytschko-Tsay  shell  elements  were  used  to  discretize  the 


side  walls  in  order  to  reduce  the  number  of  elements.  Finally,  the 
mesh  of  the  box  consisted  of  31  804  elements. 

The  Johnson-Cook  hardening  relation  [18]  was  selected  to 
model  the  aluminium  of  the  box.  There  are  more  sophisticated 
hardening  relations,  such  as  the  Zerilli-Armstrong  [49]  or  the 
Rusinek-Klepaczko  [34,35],  but  Johnson  and  Cook’s  is  probably  the 
most  widely  used  among  those  accounting  for  equivalent  plastic 
strain  ep,  equivalent  plastic  strain  rate  ep  and  temperature  ef¬ 
fects  9.  Since  numerous  efforts  have  been  made  in  the  past  to 
determine  their  parameters  for  a  large  number  of  metallic  ma¬ 
terials,  it  has  been  implemented  in  many  FE  explicit  codes.  The 
relation  is  stated  through  the  following  multiplicative  equation: 

a,  =  (,  +  B(sf)“)  (l  +  C  L»s(!))  (l  -  (£^)”)  ») 

Oo  being  the  reference  temperature  and  9m  the  melting  temper¬ 
ature.  The  required  parameters  for  the  Aluminium  6063-T5  were 
obtained  from  [19,48].  In  order  to  describe  the  material  failure  in 
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Fig.  4.  Left:  Box  model  geometry  used  for  the  analysis.  Right:  Detail  of  the  entry  wall  mesh. 


Table  1 

Parameters  used  in  the  simulation  for  the  solids. 


Material 

p  (kg/m3) 

E  (GPa) 

V 

A  (MPa) 

B  (MPa) 

n 

C 

m 

D 1 

6063-T5 

2700 

71 

0.33 

200 

144 

0.62 

0 

1 

0.2 

Steel 

7830 

207 

0.28 

- 

- 

- 

- 

- 

- 

PMMA 

1180 

3 

0.35 

- 

- 

- 

- 

- 

- 

the  zone  around  the  impact  point,  the  Johnson-Cook  model  was 
used;  the  relation  is  described  by  the  following  equation: 


£f  = 


^Di  +  D2exp^D3y- 


1  -T  D 4  In 


(6) 


where  a  and  cty  are  the  equivalent  and  the  hydrostatic  stresses, 
respectively;  the  failure  occurs  when  the  damage  parameter 
XX  A 8p/£f)  =  1.  In  this  case,  the  failure  model  was  used  with¬ 
out  taking  into  account  the  dependence  on  temperature,  strain 
rate  and  triaxiality,  and  hence  the  constants  D  2,  D 3,  D 4  and  D5 
are  set  to  zero.  Then,  failure  occurs  when  sp  reaches  the  value  D\. 

An  elastic  material  model  was  used  for  the  PMMA  window 
since  it  does  not  suffer  inelastic  deformations  nor  damage.  The 
PMMA  properties  were  obtained  from  [45].  The  projectile  was 
discretized  by  means  of  eight-node  solid  hexahedron  Lagrangian 
elements  with  reduced  integration.  The  elements  size  of  the  pro¬ 
jectile,  although  relatively  bigger  than  that  of  the  wall  elements 
in  the  impact  area,  allows  to  solve  in  a  properly  way  the  contacts 
with  the  tank  walls.  The  projectile  material  was  modelled  as  elastic 
since  no  plastic  deformation  was  observed  in  the  tests.  The  mate¬ 
rial  parameters  used  are  shown  in  Table  1. 


4.2.  Model  for  the  fluid 


It  is  expected  that,  due  to  the  HRAM  phenomenon,  the  fluid  in¬ 
side  the  tank  undergoes  too  large  deformations  to  consider  a  pure 
Lagrangian  description  as  an  appropriate  option.  For  this  reason,  a 
multi-material  ALE  formulation  with  a  second-order  accurate  ad- 
vection  method  has  been  chosen  for  the  treatment  of  the  fluid. 
Multi-material  means  that  each  element  of  the  mesh  has  the  abil¬ 
ity  to  contain  two  or  more  materials,  in  this  case  water  and  air. 
The  ALE  formulation  allows  the  motion  of  the  mesh  independently 
of  the  material  flow  without  distortion  problems. 

Fluids  inside  the  box  (water  and  air)  are  discretized  by  means 
of  eight-node  solid  hexahedron  elements  with  an  ALE  formulation 
(elform  =  11  in  LS-DYNA  notation).  Strictly,  the  fluids  are  dis¬ 
cretized  by  means  of  an  Eulerian  mesh,  but  LS-DYNA  considers  an 
Eulerian  part  as  a  special  case  of  the  ALE  formulation,  where  the 
nodes  are  fixed.  The  air  surrounding  the  box  was  also  considered, 
being  modelled  using  the  same  elements  as  in  the  fluids  inside 
the  box.  Modelling  this  air  region  is  essential  to  allow  the  water 


Slip 


mm 


Fig.  5.  Mesh  of  the  fluids. 


to  flow  into  it,  deforming  the  walls  of  the  structure.  To  achieve 
that,  water  and  air  meshes  have  to  share  the  same  nodes  at  their 
interfaces. 

Four  discretization  densities  were  analyzed  in  order  to  achieve 
an  optimal  fluid  mesh  density;  a  simplified  model  without  the  alu¬ 
minium  tube  was  used  to  perform  this  analysis.  The  results  in 
terms  of  projectile  deceleration  were  compared  with  the  classical 
solution  of  the  movement  of  a  sphere  inside  a  liquid: 

dVp  1  7 

mp  —  ~  2  PwAoCdVp  (7) 

where  mp  and  Vp  denote  the  bullet  mass  and  velocity,  pw  is 
the  fluid  density,  Aq  is  the  projected  frontal  area  of  the  projectile 
and  Q  is  a  dimensionless  drag  coefficient.  According  to  the  range 
of  velocities  considered,  a  value  of  0.4  for  Q  was  chosen  [12]. 
The  optimal  mesh  density  of  fluids,  determined  with  this  analysis, 
caused  some  problems  of  leakage  at  the  fluid/solid  interfaces.  To 
avoid  these  problems,  numerous  iterations  modifying  some  of  the 
coupling  parameters  were  made,  and  also  the  mesh  size  was  mod¬ 
ified  in  order  to  match  the  Lagrangian  one  at  the  walls  interfaces, 
obtaining  a  slightly  grosser  mesh  than  the  optimal  one  previously 
selected.  This  change  in  the  mesh  helped  to  control  the  leakage 
problems.  Finally,  the  fluids  inside  the  box  and  the  surrounding  air 
region  present  123  038  elements  (Fig.  5). 
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Table  2 

Water  and  air  parameters  used  in  the  simulation. 


A*  (5) 

Vd  (Pas) 

C(f) 

Si 

s2 

S3 

ro 

a 

c4 

c5 

Water 

1000 

0.89  •  10~3 

1448 

1.979 

0 

0 

0.11 

3.0 

- 

- 

Air 

1.22 

1.77  •  10~5 

- 

- 

- 

- 

- 

- 

0.4 

0.4 

The  water  was  modelled  using  the  following  viscous  constitu¬ 
tive  equation  ( Material  Null  in  LS-DYNA  notation) 


aij  =  2vde'i]-P8ij  (8) 

in  which  cr2J-  is  the  Cauchy  stress  tensor,  Vd  the  dynamic  viscosity, 
s-j  the  deviatoric  strain  rate,  and  the  identity  tensor.  The  devi- 
atoric  term  is  negligible  compared  to  the  hydrostatic  one,  both  for 
water  and  kerosene  due  to  the  low  value  of  their  dynamic  viscosi¬ 
ties.  The  pressure  P  is  calculated  as  a  function  of  the  compression 
pc  =  p/po  -  1,  where  p  and  po  are  the  current  and  initial  den¬ 
sities  of  the  material  respectively,  and  of  the  internal  energy  per 
unit  volume  E,  using  the  Mie-Gruneisen  Equation  of  State  based 
on  a  cubic  shock  velocity  (us)-particle  velocity  ( up ) 


P  = 


PoD2/x[l  +  (1  -  f  )/x  -  § /x2] 


Si  ^ 2 

3 2M+1 


[1-(Si  -l)/x 
for  compressed  materials  and 


S3 


JL 


+  {yo+a/i)E  (9) 


(M+i)2 


P  =  p0D2/x  +  (yo  +  a/i)E  (10) 

for  expanded  materials,  where  Si,  S2  and  S3  are  coefficients  of  the 
slope  of  the  us-up  curve.  D  is  the  intercept  of  the  us-up  curve, 
which  corresponds  to  the  adiabatic  sound  speed  on  water,  yo  is 
the  Gruneisen  gamma  and  a  is  the  first-order  volume  correction 
to  yo.  The  required  properties  and  constants  of  water  were  ob¬ 
tained  from  [8]. 

The  air  was  modelled  using  the  same  constitutive  equation 
(Eq.  (8)),  with  the  properties  of  the  air,  and  a  Linear  Polynomial 
Equation  of  State  [26].  This  equation  of  state  is  linear  in  the  in¬ 
ternal  energy  and  polynomial  in  the  compression,  and  defines  the 
pressure  P  as  follows: 


5.  Results 

In  the  following  section,  numerical  results  of  partially  filled 
tubes  (75%  and  60%)  impacted  at  velocities  of  900  and  600  m/s 
are  analyzed  and  compared  with  the  available  experimental  data. 

5.1.  HRAM  stages 

Firstly,  it  has  been  proved  that  the  ALE  model  employed  is  ca¬ 
pable  of  reproduce  qualitatively  an  HRAM  phenomenon  in  partially 
filled  tubes,  as  it  was  demonstrated  in  a  completely  filled  case  [44]. 
Each  of  the  HRAM  stages  have  been  analyzed  in  order  to  prove  that 
they  are  well  represented  in  partially  filled  tubes  simulations. 

•  Shock  phase.  When  the  projectile  penetrates  the  wall  of  the 
fluid-filled  structure,  the  impact  energy  is  transferred  to  the 
fluid  generating  a  high-pressure  hemispherical  shock  wave, 
Fig.  6(a). 

•  Drag  phase.  In  the  drag  phase,  part  of  projectile’s  kinetic  en¬ 
ergy  is  transformed  into  fluid  motion.  The  displacement  of  the 
fluid  from  the  projectile  path  generates  a  radial  pressure  field 
and  a  cavity  behind  it.  Fig.  6(b)  depicts  the  cavity  formed  in 
the  simulations. 

•  Cavitation  phase.  The  expansion  and  collapse  (oscillations)  of 
the  cavity  is  known  as  the  cavitation  stage.  In  the  model,  the 
expansion  of  the  cavity  is  well  represented,  Fig.  6(c).  The  col¬ 
lapse  is  not  captured  due  to  the  long  time  needed  for  the 
cavity  to  close. 

•  Exit  phase.  The  projectile  exits  the  tank  through  a  pre-stressed 
wall,  caused  by  the  initial  shock  stage  and  the  subsequent 
loading  by  the  fluid.  This  is  well  depicted  in  Fig.  6(d),  where 
contours  of  effective  stress  are  shown  in  the  exit  wall  before 
being  penetrated  by  the  projectile. 


P  —  C  0  +  Ci/x  +  C2/X2  +  C3/x3  +  (C4  +  C  5/x  +  C6/x2)E  (11) 


5.2.  Time  history  of  the  projectile  trajectory 


The  air  was  considered  as  an  ideal  gas  by  setting  Co  =  C\  =  C2  = 
C3  =  Co  =  0  and  C4  =  C5  =  y  -  1,  where  y  =  1.4  is  the  ratio  of 
specific  heats: 


and  the  pressure  P  is  given  by: 

p  =  (y- 1)— £  (13) 

Po 

The  properties  and  parameters  used  in  the  simulation  for  the  wa¬ 
ter  and  the  air  are  shown  in  Table  2. 

The  fluid-structure  interaction,  for  both  projectile/fluid  and 
walls/fluid,  is  achieved  by  means  of  a  penalty-based  ALE-Lagrang- 
ian  coupling  algorithm  implemented  within  LS-DYNA.  This  allows 
the  fluid  material  to  flow  around  the  structure,  but  prevents  its 
penetration  into  the  structure  mesh  applying  penalty  forces  to  the 
fluid  and  the  structure.  As  soon  as  a  fluid  node  penetrates  the  La- 
grangian  structure,  a  force  of  recall  is  exerted  to  both  the  fluid  and 
structural  node  to  make  the  fluid  node  return  to  the  surface  of 
the  structure  avoiding  penetration.  The  penalty  forces  are  propor¬ 
tional  to  the  penetration  depth  and  penalty  stiffness,  behaving  like 
a  spring  system. 


By  means  of  the  digital  high-speed  camera,  it  is  possible  to  de¬ 
termine  the  velocity  and  the  position  of  the  projectile  inside  the 
fluid.  This  experimental  data  is  compared  with  the  numerical  and 
analytical  results  [43]  (Figs.  7  and  8),  these  last  obtained  from 
Eq.  (7). 

The  figures  show  a  good  correlation  between  numerical  and  an¬ 
alytical  curves,  especially  in  those  predicting  the  position  of  the 
projectile.  The  ALE  simulations  predict  a  higher  deceleration  of  the 
projectile,  due  to  the  fluid  mesh  that  has  been  chosen  to  avoid 
the  leakage  problem.  In  general,  it  has  been  observed  that  the  re¬ 
sults  by  using  ALE  algorithm  are  very  sensitive  due  to  the  Euler 
mesh  [21,33].  Finally,  it  could  be  said  that  the  numerical  curves 
correlate  quite  well  with  the  experimental  ones  showing  the  same 
trend. 

5.3.  Cavity  evolution 

An  interesting  parameter  to  analyze,  that  has  been  less  studied, 
is  the  evolution  of  the  cavity  inside  the  fluid.  Figs.  9  and  10  show 
a  comparison  of  the  cavity  profile  obtained  both  experimentally 
and  numerically  for  different  fluid-filled  levels  on  three  different 
instants  at  an  impact  velocity  of  600  m/s.  It  can  be  seen  that  the 
size  of  cavity  predicted  by  the  simulations  is  very  similar  to  the 
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(C)  (d) 

Fig.  6.  HRAM  phases  in  a  tube  impacted  at  900  m/s  and  filled  75%:  (a)  Shock  pressure  at  t  =  0.03  ms;  (b)  cavity  formed  in  the  drag  phase  at  t  =  0.14  ms;  (c)  expansion  of 
the  cavity  at  t  =  1.2  ms;  (d)  contours  of  effective  stress  in  the  exit  wall  before  being  impacted  at  t  —  0.25  ms. 
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Fig.  7.  Comparison  of  position  of  the  projectile  vs.  time. 


size  observed  in  the  experimental  images;  this  is  of  great  impor¬ 
tance  since  the  cavity  is  the  main  cause  for  the  deformation  of 
the  walls.  Special  attention  deserves  the  portion  of  fluid  which  is 
raised  over  the  initial  fluid  level.  This  layer  is  the  main  responsible 
for  the  deformation  of  the  upper  wall. 


5.4.  Pressure  time  history 


As  was  already  mentioned  before,  two  pressure  gauges  were 
used  to  record  the  time  history  of  the  HRAM  pressure  wave  as  it 
is  propagated  through  the  fluid.  One  of  the  pressure  gauges  (PTn) 
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Fig.  8.  Comparison  of  velocity  decay  vs.  time. 


was  located  near  the  impact  point,  at  30  mm  from  the  wall  and 
75  mm  from  the  shot  line  (Fig.  3),  while  the  other  (PTf)  was 
in  the  middle  of  the  tube,  at  150  mm  from  the  projectile  tra¬ 
jectory  (Fig.  3).  The  pressure  time  histories  depicted  in  Figs.  11 
and  12  show  experimental  and  numerical  curves  representative  of 
the  pressure  measurements  at  PTn  and  PTf  at  different  velocities  in 
a  case  of  a  tube  filled  at  75%.  The  numerical  curves  correspond  to 
the  averaged  value  of  several  elements  around  each  pressure  gauge 
position.  It  can  be  seen  that  the  pressure  time  history  changes  as 
a  function  of  the  location  of  the  pressure  gauge  as  well  as  the  ve¬ 
locity  of  the  projectile.  It  is  observed  that  the  maximum  pressure 
value  is  not  accurately  predicted,  specially  on  the  PTf,  neverthe¬ 
less  the  trend  of  the  influence  of  the  different  factors,  position  and 
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Fig.  9.  Cavity  evolution  at  0.26,  0.46  and  0.60  ms  obtained  from  experiments  (top)  and  ALE  simulation  (bottom)  for  a  tube  impacted  at  600  m/s,  filled  75%. 


Fig.  10.  Cavity  evolution  at  0.13,  0.24  and  0.38  ms  obtained  from  experiments  (top)  and  ALE  simulation  (bottom)  for  a  tube  impacted  at  600  m/s,  filled  60%. 


velocity,  are  well  captured  by  the  simulations.  The  beginning  of 
the  pressure  pulse  in  the  simulations  coincides  quite  well  with 
the  experimental  data,  so  the  propagation  of  the  pressure  wave 
generated  on  the  first  stage  of  HRAM  is  well  reproduced  by  the 
simulations. 

The  impulse,  defined  as  the  area  Ap  below  the  time-pressure 
curve,  was  used  by  Wierzbicki  and  Moussa  [47]  to  characterize  the 
load  applied  by  the  fluid  in  an  aircraft  fuel  tank  under  HRAM,  un¬ 
coupling  the  fluid-structure  interaction  problem.  Fig.  13  shows  the 
impulse  at  each  pressure  transducer  as  a  function  of  volume  per¬ 
centage  for  the  different  velocities  considered.  Data  for  the  com¬ 
pletely  filled  case  has  been  added  to  the  figure  in  order  to  have 
the  complete  information.  It  can  be  seen  how  the  Ap  values  ob¬ 
tained  numerically  and  experimentally  match  very  well.  This  is  of 
great  importance  because  the  impulse,  transmitted  to  the  walls  by 
the  fluid,  which  has  influence  on  the  tube  deformation,  is  well 
predicted.  Furthermore,  this  figure  confirms  the  capability  of  the 
numerical  model  for  reproducing  the  influence  of  different  param¬ 
eters,  as  impact  velocity,  volume  percentage  and  transducer  posi¬ 
tion,  on  the  fluid  pressure. 

5.5.  Deformations 


75%,  V=900  m/s 


In  this  section,  experimental  deformation  data  in  the  entry  and 
exit  walls  of  the  tubes  are  compared  with  the  simulations  results. 


Fig.  11.  Pressure  time  history  near  (PTn)  and  far  (PTf)  from  the  impact  point  in  a 
tube  75%  filled,  impacted  at  900  m/s. 
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75%,  V=600  m/s 


Fig.  12.  Pressure  time  history  near  (PTn)  and  far  (PTf)  from  the  impact  point  in  a 
tube  75%  filled,  impacted  at  600  m/s. 


Deformation  of  the  walls  is  of  interest  due  to  its  potential  influence 
on  the  behaviour  of  contiguous  cells  that  can  be  part  of  a  whole 
fuel  tank  system.  Figs.  14  and  15  depict  the  deformation  of  the 
entry  and  exit  walls,  all  along  the  tube,  for  filling  percentage  of 
75%,  impacted  at  900  and  600  m/s.  The  time  in  the  simulations 
was  long  enough  to  assure  the  stabilization  of  the  displacements 
in  the  walls. 

The  figures  show  that  the  affected  zone  in  which  the  deforma¬ 
tion  is  produced  is  well  defined  in  the  simulations,  allowing  to 
know  the  limits  of  its  extension.  It  is  also  observed  that  the  defor¬ 
mation  obtained  in  the  simulations  is  smaller  than  that  measured 
experimentally,  nevertheless  the  trends  of  the  walls  response  to 
the  HRAM  phenomenon  are  the  same.  The  results  of  the  deforma¬ 
tion  of  the  entry  and  exit  walls  for  a  tube  filled  at  60%,  impacted 
at  900  and  600  m/s  showed  the  same  behaviour.  Considering  the 
mentioned  observations,  it  can  be  said  that  the  numerical  model 
predicts  in  an  appropriate  way  the  influence  of  the  considered 
factors,  impact  velocity  and  volume  percentage,  on  the  behaviour 
of  partially  filled  tube  walls.  This  is  important  to  identify  possi¬ 
ble  solutions  in  order  to  mitigate  the  influence  of  different  fac¬ 
tors. 


75%,  V=900  m/s 


Fig.  14.  Experimental  and  numerical  deformation  on  the  entry  and  exit  wall  of  a 
tube  filled  75%,  impacted  at  900  m/s. 


75%,  V=600  m/s 


Distance  [cm] 


Fig.  15.  Experimental  and  numerical  deformation  on  the  entry  and  exit  wall  of  a 
tube  filled  75%,  impacted  at  600  m/s. 

6.  Conclusions 

The  capability  of  the  ALE  simulation  approach  in  simulating 
the  Hydrodynamic  Ram  phenomenon  in  partially  filled  tubes  was 
studied  in  this  work.  Results  of  the  numerical  model  were  com¬ 
pared  with  experimental  results  in  partially  water-filled  aluminium 
square  tubes,  impacted  at  different  velocities. 

The  ALE  technique  is  capable  of  faithfully  reproducing  the 
stages  of  the  Hydrodynamic  Ram  in  partially  filled  tubes:  shock, 
drag,  cavitation  and  exit  phases  from  a  qualitative  and  quantita¬ 
tive  perspective.  The  cavity  evolution,  main  cause  of  the  tank  final 
deformation,  is  accurately  reproduced  as  compared  with  the  im¬ 
ages  taken  from  the  experimental  tests.  The  projectile  deceleration 
is  also  well  reproduced  for  the  different  filling  percentages;  this 
fact  is  of  great  importance,  since  the  kinetic  energy  loosed  by  the 
projectile  is  transferred  to  the  structure  (in  the  form  of  plastic  de¬ 
formation  and  damage)  through  the  fluid. 

The  pressure,  impulse  and  deformations  obtained  numerically 
have  been  compared  with  the  experimental  results  showing  that 
globally,  the  numerical  model  developed  captures  the  influence  of 
the  different  factors  considered.  This  fact  will  allow  to  analyze  in 
a  more  closely  way  the  Hydrodynamic  Ram  phenomenon  in  order 
to  try  its  attenuation  in  future  designs,  avoiding  the  expensive  cost 
of  full  scale  impact  tests. 
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